Quantum Stochastic Absorption 
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We report a detailed and systematic study of wave propagation through a stochastic absorbing 
random medium. Stochastic absorption is modeled by introducing an attenuation constant per unit 
length a in the free propagation region of the one-dimensional disordered chain of delta function 
scatterers. The average value of the logarithm of transmission coefficient decreases linearly with 
the length of the sample. The localization length is given by £ = (, m £.a/((iw + where £ w and 
£ a are the localization lengths in the presence of only disorder and of only absorption respectively. 
Absorption does not introduce any additional reflection in the limit of large a, i.e., reflection shows 
a monotonic decrease with a and tends to zero in the limit of a — » oo, in contrast to the behavior 
observed in case of coherent absorption. The stationary distribution of reflection coefficient agrees 
well with the analytical results obtained within random phase approximation (RPA) in a larger 
parameter space. We also emphasize the major differences between the results of stochastic and 
coherent absorption. 

PACS Numbers: 42.25.Bs, 71.55.Jv, 72.15.Rn, 05.40.-a 



I. INTRODUCTION 

Wave propagation in an active random medium has at- 
tracted much attention during the past decade. Recently, 
many experiments have reported lasing action of light in 
optically active strongly scattering media Q . These sys- 
tems exhibit interesting physical properties due to the 
combined effects of static disorder-induced multiple scat- 
tering and of coherent amplification/absorption [p|-^5[. 
In the extensively studied case of an electron motion 
in a random medium it is well established that quan- 
tum interference effects arising from a serial disorder in 
one-dimensional systems lead to Anderson localization 
p6| , p7t . An essential difference between electron and light 
propagation is the absence of conservation law for pho- 
tons. Light can be absorbed or amplified retaining the 
phase coherence. In most of the theoretical studies am- 
plification or absorption is modeled phenomcnologically 
by introducing an imaginary potential (optical potential) 
in the Hamiltonian. In the case of light (electro-magnetic 
waves) this corresponds to a medium with a complex di- 
electric constant. Several interesting effects have been 
predicted which include statistics of super-reflection and 
transmission [p|-[l"3[ and the dual symmetry between ab- 
sorption and amplification |lj,|l^|. Media thus modeled 
are referred to as coherently absorbing or amplifying. 

In the case of electron transport, inelastic scattering 



(due to phonons) leads to loss of phase memory of the 
wave function. Thus the motion of electrons becomes 
phase incoherent and sample to sample fluctuations be- 
come self-averaging in the high temperature limit leading 
to a classical behavior. There has been much interest in 
the effect of inelastic scattering on the coherent tunnel- 
ing through potential barriers. To allow for the possi- 
bility of inelastic decay on the otherwise coherent tun- 
neling through potential barriers, several studies invoke 
absorption | pljjl9 |. To study the above phenomenon, one 
resorts to the optical potential models (coherent absorp- 
tion models). 

In the optical potential model the potential is made 
complex V(x) = V r (x) — iVi. The Hamiltonian becomes 
non-Hcrmitian resulting in absorption or amplification of 
probability current depending on the sign of Vi . The pres- 
ence of imaginary potential (absorption/amplification) 
leads to many counter- intuitive features. In the scat- 
tering case, in the vicinity of the absorber, the particle 
experiences a mismatch in the potential (being complex) 
and therefore it tries to avoid this region by enhanced 
back reflection. Imaginary potential plays a dual role 
of an absorber and a reflector f|,|2(|. In other words, in 
such models absorption without reflection is not possible. 
Naively one expects the absorption to increase monotoni- 
cally as a function of Vi. However, the observed behavior 
is non- monotonic Pffifl . At first absorption increases and 
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after exhibiting a maximum decreases to zero as Vi — ► oo. 
The absorber, in this limit acts as a perfect reflector. 
During each scattering event an electron picks up an ad- 
ditional scattering phase shift due to Vi which along with 
multiple interference leads to additional coherence or res- 
onances in the system |2]}| . Thus, due to the presence of 
imaginary potentials, we have additional reflection and 
resonances in the system. In the presence of coherent 
absorption and quenched disorder, the stationary distri- 
bution for reflection coefficient has been calculated 0. 
This has been done within random phase approximation 
(RPA) using the invariant imbedding method p2| . The 
stationary distribution is given by: 

\D\exp{\D\)exp(-^) f < . m 
P s (r) = _ ^ 2 for r < 1 (1) 

= 



for r > 1. 



Here D is proportional to Vi/W, W being the strength of 
disorder. Notice that the distribution has a single peak 
which shifts towards r = with increasing absorption 
strength Vi. However, the exact distribution obtained 
numerically for strong disorder and strong absorption 
shows significant qua litative departure from this analyti- 
cal distribution ||,jl3) . For sufficiently strong absorption, 
the numerically obtained stationary distribution shows a 
double peak structure. In the limit Vi — > oo the distribu- 
tion becomes a delta function at r = 1. This corresponds 
to the limit where the absorber acts as a perfect reflector. 

To this end we would like to develop a model where 
absorption does not lead to concomitant reflection and 
additional resonances as discussed above. Recently, such 
a stochastic absorption model was developed by Prad- 
han based on the work of Biittiker |^5|,^6| . In his 

treatment several absorptive side-channels are added to 
the purely elastic channels of interest. A particle that 
once enters the absorbing or the side-channel never re- 
turns back and is physically lost. He has obtained the 
Langevin equation for the reflection amplitude R(L) for 
a random medium of length L by enlarging the S'-matrix 
to include side-channels. In continuum limit the equation 
for R(L) is [|H|: 



dR 
dL 



— = -aR(L) + 2ikR(L) + ikV(L)[l + R(L)} 2 , (2) 



where a is the absorption parameter and V(L) is the 
random potential representing the static disorder. Inter- 
estingly, within the random phase approximation (RPA) , 
the stationary probability distribution for the reflection 
coefficient P s (r) ( for L — ► oo ) is again given [ p3|j24| 
by Eq.0. In our present work we develop another simple 
model for absorption which can be readily used to study 
the case of amplifying medium as well. The medium 
comprises of random strength delta function scatterers 
at regular spatial intervals a. To model absorption (leak- 
ing out) of electrons, an attenuation constant per unit 
length a is introduced. Every time the electron traverses 



the free region between the delta scatterers, we insert 
a factor exp(-aa) in the free propagator following Ref. 
p7[ . We find that this method of modeling absorption 
does not lead to additional reflection and resonances as 
in the case of optical potential models. We obtain the 
localization length and study the statistics of reflection 
and transmission coefficients. The stationary distribu- 
tion of reflection coefficient agrees with Eq.llin a larger 
parameter space. Following earlier method [23], the con- 
tinuum limit of our model leads to the same Langevin 
equation (Eq. ||) for R(L) where a is replaced by 2a. 
Naturally, agreement of our result with Eqjl] follows. In 
Sec.[n]we give the details of our model and the numerical 
procedure. The section after that is devoted to results 
and discussion. 



II. THE MODEL 

We carry out calculations on the wave propagation 
in an absorbing medium characterized by an attenua- 
tion constant a and interspersed by a chain of uniformly 
spaced independent delta-function scatterers of random 
strengths. The i th delta-function scattering center is de- 
scribed by a transfer matrix |2§f| 



M; 



- iqi/2k 
iqi/2k 



-iqi/2k \ 
f iqi/2k ) 



where qi is the strength of the i th delta-function. The 
qi's are uniformly distributed over the range — W/2 < 
qi < W/2, i.e., P(q t ) = l/W. Here W is the disorder 
strength. We set units of h and 2m to be unity. The 
energy of the incident wave is E = k 2 . For further anal- 
ysis, W and a are scaled with respect to a and are made 
dimensionless. Propagation of the wave in-between two 
consecutive delta- function scatterers separated by a unit 
spacing (a = 1) can be described by the matrix 



X 



^ik—ct 

e~ lfcH 



The total transfer matrix for the L-site system is con- 
structed by repeated application of Mi and X [ pSfl : 

M = M L X....XM 2 XMi. 

From M the reflection and transmission amplitudes are 
calculated using 



and 



R 



T = — 



M(2, 1) 
M(2, 2) 



detM 

M(2,2)' 



2 



The reflection and transmission coefficients are r = \R\ 2 
and t = \T\ 2 respectively and the absorption is given by 
er = 1 — r — t. Thus, due to absorption the total flux is 
not conserved and we have r + 1 =/= 1 . 



III. RESULTS AND DISCUSSION 

In our studies we consider at least 10,000 realizations 
for calculating various distributions and averages. In the 
case of stationary distributions, the length of the sam- 
ples considered were about 5 to 10 times the localization 
length. We also verified that the corresponding distribu- 
tions or averages do not evolve any further with increas- 
ing sample length L. All results are shown for incident 
energy E = k 2 = 1.0 unless specified otherwise. 

We first consider the behavior of (Int). The angu- 
lar bracket denotes the ensemble average. In Figjl] we 
plot (Int) as a function of length L for ordered absorp- 
tive medium ( W = 0.0, a — 0.05 ), ensemble av- 
eraged disordered non-absorptive medium ( W = 1.0, 
a = 0.0 ) and disordered absorptive medium ( W = 1.0, 
a = 0.05,0.1,0.15 ). In all the cases transmission decays 
exponentially with the length. The absorption-induced 
length scale £ in random medium associated with the de- 
cay of transmission coefficient is always less than both 
£ a and £ w . The localization length for disordered nonab- 
sorptive medium scales as ^9) £ w — 96k 2 /W 2 and for or- 
dered absorptive medium, as £ a = 1/a. In Fig. |] we show 
the plot of l/£ versus l/£ w + l/£a obtained by changing 
a for various values of disorder strength W. We have nu- 
merically calculated l/£ for the different cases. All the 
points fall on a straight line with unit slope indicating 
the relation l/£ = l/£ w + l/£ a - Such a relation ex- 
ists for the case of coherently absorbing and amplifying 
media §,fL§. 

To study the nature of fluctuations in the transmission 
coefficient, in Fig.|| we plot, on log-scale, average t ({*)), 

root-mean-squared variance t v = \J (t 2 ) — (t) 2 and root- 
mean-squared relative variance t rv =t v j (t) as a function 
of length L for W = 1.0 and a — 0.01. We see from the 
figure that t v is less than (t) and t rv is less than unity 
upto L/£ w 3. But, beyond that t v becomes greater than 
(t) and t rv crosses unity making the transmission coef- 
ficient a non-self-averaging quantity. This implies that 
the fluctuations in t over the ensemble of macroscopically 
identical samples dominates the ensemble average. The 
transmission coefficient becomes very sensitive to spatial 
realizations of the impurity configurations. 

In Fig. |] we show the distribution P(t) at different sam- 
ple lengths for a = 0.01 and W = 1.0. For small lengths 
L, resonant transmission dominates and P(t) peaks at 
a large value of t. In fact for L — > 0, P(t) — » S(t — 1). 
As the length becomes comparable to the localization 
length L ~ £, multiple reflections start dominating. Con- 
sequently, the time spent inside the medium increases 
leading to more absorption. Thus, the peak of the distri- 



bution shifts to smaller values of t and the distribution 
broadens due to randomization by disorder. In the long 
length limit L » £, the distribution develops long tails 
and its peak shifts towards t = 0. The transmittancc 
shows large sample-to-sample fluctuations and becomes 
a non-self-averaging quantity. It is the tail of the dis- 
tribution that determines the behavior of fluctuations. 
Finally, as expected, for L — > oo , P(t) — > 8(t). 

From the previous discussion it is clear that the trans- 
mittance becomes non-self-averaging due to the appear- 
ance of tails in the ^-distribution. These tails owe their 
existence to the presence of resonant realizations in the 
ensemble. Hence, it is worthwhile to investigate the na- 
ture of resonances and the effect of absorption on them. 
Specifically, we would like to understand if the presence 
of absorption would give rise to any new resonances. It 
is well known from the studies in passive disordered me- 
dia that the ensemble fluctuation and the fluctuations 
for a given sample as a function of chemical potential or 
energy are expected to be related by some sort of ergod- 
icity ]30|] , i.e., the measured fluctuations as a function 
of the control parameter are identical to the fluctuations 
observable by changing the impurity configurations. In 
Fig.^(a) we show the plot of t versus k for W = 1.0 and 
a = at L = 100 for a given realization of the random 
potential. Figure ||(b) shows a plot of t versus k for the 
same realization but with a = 0.01. By mere visual in- 
spection one can see that the only effect of absorption, 
apart from reducing the value of transmission, is to in- 
crease the width of resonance peaks for the passive case. 
Thus the presence of absorption does not introduce any 
new resonances. This can be seen from Fig.||(c) and (d) 
which emphasize FigJ^(a) and (b) respectively by enlarg- 
ing a narrow region between k = 1.5 to k = 2.0. Wc 
do not see any new peaks in the transmission spectrum 
for absorptive case. Similar effect is observed in case of 
reflection also. 

We now turn our attention to the statistics of reflec- 
tion coefficient. In Fig|j| we plot (Inr) as a function of 
length L for a fixed value of disorder strength W — 1.0 
and different values of absorption strength a as indicated 
in the figure. It increases with L initially and for L>^, 
it saturates. At any L, (Inr) \w,a < (Inr) \w,o- This is 
in contrast to the behavior observed for the case of co- 
herent absorption. As we know, in the case of coherent 
absorption, the reflection coefficient tends to unity for 
absorption strength becoming very large. In this regime, 
the predominantly reflecting nature of optical potential 
makes reflection larger than that in the corresponding 
passive case. 

In Fig.^ we have shown P s (r) for various values of a. 
In the small a range, i.e., for a = 0.001, the distribu- 
tion has a peak at large r. As we increase a the peak 
shifts to smaller values of r. The thick line shows the fit 
obtained using the analytical expression given in Eqjl]. 
In the limit of large a, the distribution tends to become 
a delta function at r = 0. This is in sharp contrast to 
the behavior observed for coherent absorption 0. The 
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distribution is always single peaked. For all non-zero val- 
ues of a the medium acts as an absorber only and there 
is no additional reflection due to absorption. Figure ^ 
shows a monotonic decrease of saturated value of aver- 
age reflection coefficient (r) as a function of a. The 
average absorption, defined as (a) = 1 — (r) — (t), in- 
creases monotonically with increasing a and in the limit 
of a — y oo saturates to unity in contrast to the optical 
model wherein it tends to zero. In the case of the opti- 
cal model, the absorption coefficient is a non-monotonic 
function of absorption strength U and for values of U 
near the peak the stationary distribution of reflection co- 
efficient displays a double peak PJH| . In fact our model 
exhibits the properties in agreement with physical expec- 
tations of an absorbing medium, i.e., stronger the absorp- 
tion lesser are the reflection and transmission across the 
medium. 

Finally, we discuss the phase distribution. Figure ^ 
shows the stationary distribution of phase of the reflected 
wave for a fixed disorder strength W = 1.0 and various 
for values of a. For small values of disorder one generally 
expects the phase distribution to be uniform if the sys- 
tem size is around the localization length. This is seen in 
Fig.|(a) for the case of weak absorption. As we increase 
a the phase distribution develops two distinct peaks - 
a feature observed for coherent absorption also. This is 
related to the fact that the localization length decreases 
with a. We would like to point out that the stationary 
distribution P s {r) is same within the RPA for the case 
of coherently as well as stochastically absorbing media. 
Unlike in the case of coherently absorbing medium, Eq.[l] 
seems to be valid in a larger parameter space for the 
stochastic absorbing medium where the RPA may not 
be valid. The parameters for validity of the RPA are 
determined by the observation of uniform phase distri- 
bution. However, beyond the RPA, P s (r) for the case of 
stochastic and coherent absorbing media are qualitatively 
distinct from each other. 

In conclusion, we have studied a new model of quan- 
tum stochastic absorption. The behavior observed for 
transmission and reflection coefficients is in accordance 
with physical expectations of an absorbing medium. This 
model can be extended to the case of stochastically am- 
plifying medium. It exhibits duality between absorption 
and amplification which has received much attention re- 
cently. Results for this will be reported elsewhere |Uj . 
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FIG. 3. Average transmission coefficient (t), variance of t 
(t v ) and relative variance of t (t rv ) as a function of L/£ for 
fixed disorder W = 1.0 and absorption a = 0.01. 
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FIG. 1. Average of logarithm of transmission coefficient t 
versus length L for different values of absorption strength a 
and disorder strength W as indicated in the figure. 
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FIG. 4. Distribution of transmission coefficient t from a 
disordered absorptive system with W = 1.0 and a = 0.01 at 
different lengths L. 




FIG. 2. l/£ versus 1/£ TO + l/£a, where ^ w is the localiza- 
tion length for non-absorptive disordered system, £ a is the 
localization length for ordered absorptive system and £ is the 
localization length for absorptive disordered system. 
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FIG. 5. Transmission coefficient as a function of incident 
wavenumber k for (a),(c) disordered non-absorptive sample 
W = 1.0 of length L — 100 and (b),(d) disordered absorptive 
sample W = 1.0, a = 0.01 of length L = 100. 




FIG. 6. Average of logarithm of reflection coefficient ver- 
sus sample length for a fixed value of disorder W = 1.0 and 
different values of absorption a indicated in the figure. 
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FIG. 7. Stationary distribution of reflection coefficient for 
a fixed value of disorder strength W — 1.0 and different val- 
ues of absorption a. The thick line shows the single parame- 
ter fit of analytical expression Eqn.ji] with (a) D = 0.197 for 
a = 0.001,(b) D = 1.92 for a = 0.01 and (c) D = 20.53 for 
a = 0.1. 
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FIG. 8. Average value of reflection coefficient ((r)), trans- 
mission coefficient ((£)) and absorption ((<j)) versus absorp- 
tion strength a for W = 1.0 and L/£ = 10. 
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FIG. 9. Stationary distribution of phase of reflected wave 
for fixed disorder strength W = 1.0 and different values of 
absorption a. 
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